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The hypernuclear Auger effect, given by the de-excitation of a A hypernucleus by means of the 
transition of a A hyperon from an initial to a lower lying final single particle state in conjunction 
with neutron emission from the core nucleus is studied in relativistic DDRH field theory. Baryonic 
interactions are obtained from the Bonn NN potential by Dirac-Brueckner theory and theoretically 
derived scaling laws for the meson-hyperon vertices. The model is applied to the A'^Pb hypernucleus. 
A and nucleon bound states as well as scattering states are calculated self-consistently in mean-field 
approximation. The Auger spectra of the emitted neutrons are of complex structure due to a 
huge combinatorial number of possible hypernuclear transitions. The sensitivity of the neutron 
Auger spectra to changes in the Act and Aoj vertices is investigated. The theoretical results show 
that experimental applications of the hypernuclear Auger effect will require special efforts, e.g. by 
tagging on the energy of the initially created hyperon. 
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Cs| ■ I. INTRODUCTION 

> 

' In order to obtain a more profound understanding of the interactions between baryons with and without strangeness 
, the major goal of hypernuclear physics is to explore the spectralproperties of A hyperons in a nucleus. The experimen- 
■ tally well confirmed existence of A single particle states |[ ^ ^, ^, ^ indicates the importance of static mean-field 
interactions. While the overall features are reasonably well reproduced by non-relativistic |Q, || and relativistic ap- 
(N proachcs §, [u) open questions exist on finer details of hypernuclear spectra. In light nuclei, A < 16, the smallness 
' of the spin-orbit splitting and other fine structures of A single particle spectra are well described by multi-configuration 
I shell model calculations with empirical matrix elements |l2[ . These calculations emphasize the importance of residual 
interactions and configuration mixing. Theoretical studies of dynamical self-energies in A nuclei, predicting rather 
strong contributions from the coupling of the A to excitations of the nuclear core, point in a similar direction ||. 
This, however, is not fully complying with recent observations of well separated and relatively sharp A single particle 
^ ^ structures in the medium-mass ^^V and ^^Y hypernuclei ^ |^. The measured spectra are compatible with a much 
stronger spin-orbit strength than expected from the low mass region, amounting to spin-orbit interaction energies of 
1 to 2 MeV, but the data do not indicate a strong damping of the A single particle states. 

The conclusion to be drawn from these - partly conflicting - findings is that there is a clear necessity for more 
$H ' precise measurements and an enlarged body of data. For more decisive results on A-nucleus interactions medium and 
heavy mass nuclei can be expected to be much better suited. In light nuclei threshold effects from the weak binding of 
A states and finite size effects contribute strongly to the observed spin-orbit splitting, thus not giving direct access to 
the wanted information on the genuine A interaction strengths ||l^ . The calculations in rcf. jl4| indicate that these 
effects decrease rapidly beyond the silicon mass region. A standard tool for spectral investigations is 7-spectroscopy, 
detecting the photons from the hypernuclear transitions of the A particle after it was created in an excited state e.g. 
in a A{Tr, K)\A reaction. In a recent {K~ ,tt~) experiment on ^^C at BNL/AGS progress has been made by achieving 
a much better energy resolution than before [ p^ . These measurements, by the way, confirm the small A spin-orbit 
splitting in light nuclei. By technical reasons 7-spectroscopy is not applicable in medium and heavy mass nuclei. 

A promising alternative for spectroscopy in medium and heavy mass hypernuclei is the observation of Auger 
neutrons, emitted during the de-excitation of the hypernucleus after the initial creation of a A in an excited state. 
The hypernuclear Auger process was discussed first by Likar et al. |l6[ . More recently, the idea has been revived 
and worked out in much more detail for an experiment proposal at JLAB ||l7|| . The hypernuclear Auger effect is the 
direct de-excitation of an excited A single-particle state in a hypernucleus by the emission of a neutron. The process 
is analogous to the well known atomic Auger effect where the de-excitation of an electron kicks out another less bound 
one. In the case of the hypernuclear Auger effect the energy spectrum of the emitted neutrons reflects the A singe- 
particle level structure, but folded with the neutron single-particle spectrum. The mechanism of the hypernuclear 
Auger effect is schematically shown in flgure |l|. In order to occur the separation energy of at least the A Is-orbit must 
be larger than the separation energy of the valence neutron. Therefore, because of the much weaker A binding, the 
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FIG. 1: Schematic picture of the hypernuclear Auger effect: A A hyperon is de-excited into a lower single particle level,thereby 
transferring energy and momentum to a valence neutron which is emitted through this process. 



hypernuclear Auger effect appears with a strength sufficient for measurements only in intermediate and heavy mass 
nuclei where neutron and A separation energies of comparable magnitude are available. Already from this simplified 
picture it is clear that a single transition of the A will yield a wealth of peaks in the neutron's energy distribution, 
requiring a very careful analysis. For the proposed JLAB experiment JT?! a neutron energy resolution of better than 
50 keV is envisaged. As nuclei to be studied in the experiment Ph and \J are proposed. In this mass region it can 
yield valuable information about the A single particle spectra, possibly even resolving finer details as the spin-orbit 
splitting. One purpose of this paper is to investigate in a realistic model predictions for spectral distributions which 
may be used to estimate constraints on future measurements. 

In this paper we report on our calculations of the hypernuclear Auger spectra within the density dependent rela- 
tivistic hadron field theory (DDRH). In section || we describe the theoretical description of the hypernuclear Auger 
effect in the relativistic DDRH approach with density dependent meson-baryon vertices. We shortly sketch the main 
features of DDRH and the numerical realization. In section III results are presented and possibilities to extract 



information from the extremely complex spectra are discussed. In section IV we summarize, discuss our results and 
draw some conclusions. 



II. THE HYPERNUCLEAR AUGER EFFECT 



A. Formalism 



The emission of a neutron by the A de-excitation of the hypernucleus is described as a decay of an initial A single 
particle state into a 2 particle- 1 hole configuration where the A is coupled to a neutron particle- hole core excitation 
with an energy above the particle emission threshold. The nucleon and hyperon single particle states are obtained in 
relativistic mean-field (RMF) approximation. Residual interactions among the neutron particle-hole configurations 
are neglected, i.e. an uncorrelated quasiparticle description is used. 

For the transition operator V the one-boson-exchange parametrization of the full Dirac-Brueckner (DB) G-matrix 
with density dependent vertex functionals r(/5) pO| is used. It is the same interaction as applied in the structure 
calculation. Since the A hyperon is electrically neutral and an iso-scalar particle only the a and the lo mesons 
contribute: 

y = Fo-aC/Oa)^ ^'^aN^PN) - F^jaIpa)^ lon{Pn) (1) 

Taking |0) to be our many-body ground state, assumed here as the 0+ ground state of a spherical nucleus, the 
initial state is in second quantization formulation given by aj^ |0), where aj^ is the creation operator for a A state 
with the set of quantum numbers A^. In the final state the hyperon is attached to a particle- hole excited nuclear 
core. The excess energy and momentum is carried away by the emitted neutron occupying an unbound single particle 
continuum state in the nuclear mean-field potential: 
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is the particle-hole excitation operator with angular momentum J„^ , M„^ defined through: 



(3) 



where a = (— )-'+™aj__TO denotes a hole creation operator. In eq.^ the a\^^ and Aj^^ are coupled to total angular 
momentum j^, mp. 

The diff'erential widths dTj^^^ describing the decay of an initial state is determined by the transition matrix 
elements of the AA^ interaction V: 
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(4) 



where 



{71} = { 



"^A„ , TO/3 , J/3 , jAfl , Jnn , J„- 1 



_ I indicates the incoherent summation over degenerate initial and final 



sub-states, including the appropriate phase space factors |18| due to the neutron emerging with momentum k. M is 
the mass of the initial hypernucleus. The orthogonality and completeness relations of the Clebsch-Gordan coefficients 
Ip^ allow to convert this expression into the equivalent form of an incoherent sum over matrix elements of uncoupled 
states: 



{72} 



M2 
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(5) 



where {72} = |mA„ , mA^, , m„^ , "i^-i , Ja^ , Jn^ , Jr,^ 
determined by nucleonic particle-hole fluctuations of the nuclear mean-field absorbing the energy-momentum transfer 
from the A transition. This structure becomes even more obvious by expressing the matrix elements in terms of the 
appropriate non-diagonal elements of A and nucleon one-body density matrices in momentum space: 



From equations and (||) it is seen that the Auger process is 



{0\aAf,ani,Val^_^al^ |0} = / d'^q pA'A{q)V{q)pn-in{q) 
Here, the transition densities pij are given by: 
p,,{q) = [ d^xe''^-^,{x)Hjix) 



(6) 



S ((qO - [E, - E,)) I dn dr | 5](-)^n;(g)lA^(r)jA(gr) ] V'.(r, f)f V-, (r, r) (7) 

/I, A 

r is either 1 or 7'' for a- and w-exchange, respectively. Eij are the single particle energies of the states ipij. In 
the second line the spatial part of the plane wave is expanded into partial waves. A more detailed description of the 
evaluation of the matrix elements is given in appendix The wave functions and single particle energies Eij are 
taken from a DDRH calculation which is described in section II B 



B. The DDRH Model 

The density dependent relativistic hadron field theory (DDRH) ||2^, |2l[ is a relativistic Lagrangian field theory with 
baryons and mesons interacting by density dependent coupling functionals. The meson-baryon interaction part of the 
DDRH-Lagrangian used here is chosen as in |20) : 

Ant = *r,(/5)*a-^r^(p)7^*c^^ 

-iwfp7^*p^ - c^fQj^.'^fA'' (8) 

including isoscalar interactions from the scalar a and the vector uj meson and isovector vector interactions given by 
the yO-meson. The last term accounts for the electromagnetic interaction where Q is the charge operator. The density 
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dependence of the baryon-meson vertices is determined self-consistently in mean-field approximation from Dirac- 
Brueckner {DB) self-energies by first solving the Bethe-Salpeter equation for the in-medium two-particle scattering 
amplitude and then extracting the medium modified vertices by expanding the DB Hartree-Fock self-energies in terms 
of the meson-exchange self-energies obtained from equation (g) |2^. Medium modifications contribute at various 
levels, e.g. the dressing of the baryon propagators by self-energies and the Pauli-blocking of intermediate states. 
A Lorentz invariant and thermodynamically consistent field theory is retained by expressing the density dependent 
vertices in terms of functionals of Lorentz scalar combinations pi of the baryon field operators ip [ pO| . The description 
of infinite nuclear matter and pure isospin nuclei within DDRH theory has been shown to be very good using the 
Bonn A and the Groningen NN-potentials |^0[ |T] . 

In the extension of the DDRH theory to the strangeness sector and applications to hypernuclear structure 
calculations were discussed. Since DB calculations describing the in-medium interactions of the complete baryon 
octet are not available we follow here the approximation scheme discussed also in . The essential step is obtained 
from the diagrammatic analysis of DB hyperon self-energies, showing that the coupling functionals of nucleons and 
hyperons are related in leading order by a simple scaling law given by the ratio of the free space meson-nucleon and 
meson- hyperon Born-term coupling constants gNa and gva, respectively 

^Yaipy) « ^^'^NaiPr) = RYa^Naipv)- (9) 
gNa 

The scalar scaling factor Rf^a was taken from the Jiilich-model (i?Aa- — 0.49) (23), including explicitly tttt and KK 
dynamical correlations in the 0+ scalar-isoscalar meson channel. The scaling factor for the AAw vector vertex is 
obtained phenomenologically from a least-square fit to A separation energies, resulting in i?Ai^ = 0.553 |^^ . However, 
in ref. it was noticed that the binding energies in the light nuclei below oxygen seem to behave differently from 
the systematics obtained from A > AQ. The A separation energies in the light nuclei are described with the same 
accuracy only if the repulsion from the lo meson is reduced by about 5% to Rf^^ = 0.542. DDRH results for single A 
hypernuclei are compared to the presently available set of world data in . An agreement on the percent level is 
obtained, being at least of the same quality as those of purely phenomenological mean-field models. 



III. RESULTS 



In this work we consider the hypernucleus ^''^Pb as representative for the heavy mass hypernuclei (it is one of 
the nuclei that are going to be studied in the JLAB experiment). The wave functions and single particle energies 
used for the evaluation of the matrix elements are calculated self-consistently by solving the DDRH field equations in 
relativistic mean-field (RMF) approximation. Details of the numerical approach and the model parameters are given 
in [ p^ . For the present application we need information on knocked-out, unbound neutron states. For that purpose 
the single particle continuum was discretized by enclosing the system in a huge box of size R = 150/m. Since the 
spacing of the discretrized continuum levels behaves as ^ 0{1/R^) the use of such a large quantization volume ensures 
quasi-continuous energy spectra for the neutron scattering states, allowing to resolve single particle resonances and 
other continuum structures resulting from the calculations. 

The continuum wave functions are calculated in the self-consistently obtained ground state mean-field potentials. 
Thus, final state interactions are taken into account on the level of static mean-field self-energies. The approach assures 
orthogonality of bound and unbound wave functions thus avoiding the unphysical non-orthogonality contributions 
inherent to phenomenological approaches. 

Since the previous applications of the DDRH theory show that the experimental A spectra are reproduced especially 
accurate in heavy nuclei our approach is well suited for the calculation of the Auger neutron spectra from ^"^^Pb. 
Numerical values of the single particle energies for occupied neutron states and bound A states are displayed in 
table I 

Because of the high level density and moderate separation energies heavy hypernuclei are most suitable for the 
Auger effect. However, at the same time these apparent advantages are, unfortunately, a potential source of problems 
for experimental work. The huge amount of combinatorial possibilities for transitions, illustrated in table ^ and 
indicated in figure || for the case of initially populating the ly^-shell in ^"^Pb, leads to Auger spectra of a rather 
complicated shape making in many cases an unambiguous identification of transitions and assignment of quantum 
numbers almost impossible. 

The problem is apparent from figure ^ where the complete neutron emission spectrum from the ^^Pb hypernucleus 
is displayed, summed over all energetically open A levels. Experimentally, spectra of a similar structure have to 
be expected. In addition, state-dependent weighting factors from the production vertex of the initial A state will 
be superimposed. From figure |^ and figure ^ it is obvious that experiments will be confronted with spectra of high 
complexity. Before observables of physical interest can be accessed the data will have to be analyzed in a more selective 
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TABLE L DDRH results for A and neutron single particle energies in a Pb entering into the Auger-calculations. The standard 
set of interaction parameters was used |l4|. 
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FIG. 2: Level schemes of the A hyperons and neutrons that 
are involved in the hypernuclear Auger effect. The levels are 

displayed in the physical scale as shown in table |. The bars piQ. 3: Full spectrum of the Auger transition strengths in 

mark the single particle levels involved by the de-excitation of 209 including contributions by initial population of the Ih, 

an initial Ig-shell A. The maximum energy, released when the ig^ ^s, 2d, 2p, 2s and If A orbitals and their subsequent 

A drops down to the Is-orbit, allows to emit neutrons from de-excitation by neutron emission, 
the lii3/2 valence orbit down to the 2^3/2 or 3si/2 shells for 
the 1(77/2 and lgg/2 initial states, respectively. 



approach. From the discussion it is obvious that much of the structure will be produced by incoherent superpositions 
of contributions from the variety of orbitals in which the A was inititally produced. By a precise energy tagging of the 
produced kaon and the outgoing electron (for the case of electromagnetic production of hypernuclei, as will be the case 
at JLAB) it might be possible to determine the initial A-state accurately enough and that only the Auger-neutrons 
related to the de-excitation of that specific state can be recorded selectively in a coincidence measurement. 

In this context, it is of interest to consider the Auger transitions on the limits of the energetically accessible range 
of states. Assuming energetically sharp states, i.e. the A and neutrons are in good quasiparticle configurations and 
damping effects are negligible, the present calculation predicts that the Auger process can only take place if the A 
initially is produced in the 2s'^-orbit or above (see table |), because otherwise the A transition energies are less than 
the lowest neutron separation energy. Next to the 2s^-orbit we find the doublet of 1/^ states which is of interest 
because it allows to observe the A spin-orbit splitting, at least in principle. Since the energy window available from 
populating the 1/^-orbits is still rather narrow the resulting Auger neutron spectrum is of a comparatively simple 
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TABLE II: Transitions contributing to the de-excitation of the Ig'^-states. Transition widths for the A initial states lfl9/2 and 
lg7/2 are denoted by Tigg^^ and Fig^^^, respectively. 
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FIG. 4: Decay widths of the Auger transition strengths in a Pb with an initial If A state. 
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structure. In figure ^ results for the Auger spectra produced by the I/7/2 and I/5/2 states are compared. It is seen 
that in both cases only a small number of final states occurs. Even more, only the transitions to the lji3/2 neutron- 
hole final state yield a significant strength. Thus it is even possible to clearly obtain the spin-orbit splitting of the 
1/^-shell. The prominent strength of this doublet is also fairly model independent. Considering the kinematically 
allowed phase space of this transition only, one might expect that wave function effects which are sensitive to the 
details of the interaction could strongly influence the relative strengths between the emission of the li- or 3p-shell 
neutrons, which are in energy almost degenerate (see table ^. Nevertheless, due to the high degeneracy of the 1*13/2 
neutron-orbitals the multiplicity of these neutrons will be greatly enhanced. Therefore one can expect to observe in 
the de-excitation spectrum of a 1/^-state a clear line doublet belonging to the two spin-orbit partners of the If^ shell 
falling down to the Is orbital and thereby knocking out the lii3/2 neutron. 
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-20- 
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2p- 



lf= 

2s- 
ld= 



Ip- 

Is- 



C=0.6 



C=1.0 



C=1.3 



FIG. 5: A single particle energy spectra in j^^Ph for the three different sets of coupling constants discussed in the text. 



The general case is, however, more ambiguous. The doublet structure which one might expect as the signature 
of the A spin-orbit splitting does not always show up in the spectral distribution since in many cases the matrix 
elements depend sensitively on binding energies and other wave function effects. As a representative case, we study 
the influence of the spin-orbit interaction strength on the Auger neutron spectra in more detail for the Ig^ shell. This 
case is well suited, since the spin-orbit splitting is sufficiently large and the spectrum offers already some complexity. 
The possible transitions for this configuration together with the corresponding transition rates are shown in table ||. 

The dependence of the Ig^-shell spectra on variations of the spin-orbit strength of the A-nucleus potential is 
investigated by changing the relative and the absolute coupling strength of the a and the uj meson to the A, keeping 
the overall single particle structure of the A spectrum fixed. Numerically, this is realized by observing that in our 
relativistic mean-field theory with scalar and vector self-energies and Ucj, respectively, the leading order non- 
relativistic Schroedinger-type central potential is given by Uo — U^i — Ua and the strength of the spin-orbit potential 
is determined by Uis = Ui^ + U^. Hence, we can relate spectral effects from variations of the spin-orbit strength by a 
factor C to a scaling of the scalar and vector self-energies Ua and , respectively, according to 

DUo = pU^-aU, . (10) 

We preserve the self-consistency between the nuclear mean-field and the underlying baryon-baryon interactions by 
scaling the Act and Auj vertices, equation (^, by the same factors a and respectively. 

Choosing C as our external parameter and using always D = 1 as a constraint, the spin-orbit splittings of the A levels 
can be varied over wide ranges while keeping the changes in the overall structure of the A single particle spectrum on 
a minimal level. The A single particle spectra obtained for C = 0.6, 1.3 (corresponding to {a, (3) = (0.6, 0.5), (1.3, 1, 4), 
respectively,) are compared to the results for C = 1 in figure |. For C = 0.6 the spin-orbit splittings are reduced by 
about a factor of 2 and an increase of about the same size is found for C = 1.3. These (strong) variations will surely 
cover the full range of uncertainties about the A spin-orbit potential. 
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FIG. 6: Comparison of the spectral Auger strength distributions produced by transitions from the Ig^ shell for interactions 
with different spin-orbit splittings. The Auger transition widths F are shown as functions of energy En of the outgoing neutron. 
The upper and lower arrows indicate transitions from the lgg/2,7/2 to the IP3/2 and lpi/2 final states, respectively. 



The spectral distributions of Auger transitions widths F, obtained by putting the initial A in 139/2,7/2 orbits are 
displayed in figure ^. Compared to the standard case C = 1 the spectral structures and transition strengths are 
changed drastically when using C — 0.6 and C = 1.3. As an overall feature the calculations show a concentration 
of strength in a few states for C = 1 and C — 1.3 while a more equilibrated distribution is found for C = 0.6. 
The apparent pile-up of strength at low neutron energies for the normal and strong coupling cases are related to 
the variations in the values of overlap matrix elements due to changes in binding energies and single particle wave 
functions. Analyzing the dynamical content of the response functions by calculating sum rules for various moments 
of the excitation energy one finds a disappointing small sensitivity on the spin-orbit interaction strength. Hence, it is 
unlikely that spin-orbit effects will contribute significantly and on an observable level to the total spectral strength. 

On the level of individual transitions there are, however, signals for spin-orbit effects visible. In figure O the 
transitions Ig^^j 7/2 ~* ^P3/2 indicated by the upper arrows while lower arrows denote Ig^^j 7/2 ~^ 
length of the arrows corresponds to the 15^/2 ^ I57/2 ^^^^SY splitting which obviously depends directly on the spin- 
orbit interaction strength. In addition, details of the spectral distributions, e.g. the clustering of strength in certain 
energy regions, also depends on the overlap of wave functions by which the transition matrix elements and therefore 
the transition widths F are determined (see equations (||) and (^). 

It might seem somewhat disadvantageous to start out doing Auger experiments with a system of the complexity 
of a heavy nucleus like lead. However, aiming at resolving the A spin-orbit structure in heavy hypernuclei, which is 
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still a controversial question as was mentioned in the introduction, one needs to find hypernuclei that provide a core 
nucleus of very low spin. In the ideal case this should be J'^ — 0+ in order to eliminate - or at least to suppress - 
effects from A angular momentum - core nucleus spin interactions leading to an additional splitting superimposed 
on the pure A spin-orbit splitting. In this respect the frequently used ^^Y and ^^V hypernuclei are very unfavorable 
cases because the ^^Y and core nuclei carry ground state spins of 4~ and 6+, respectively! The broadening of the 
1-orbits observed in these nuclei [1] does not exhibit a proper I ■ s systematics which might be related to contribution 
from residual interactions involving the huge core spins [25] . Hyperisotopes matching the requirements of having a 
0+ core ground state and appearing with at least around 10% natural abundancy would be ^^Zr and ]^^^Sn. Lighter 
isotopes will gradually cease to exhibit Auger de-excitations. 

IV. SUMMARY AND CONCLUSIONS 

The Auger-neutron transition rates for the de-excitation of single A hypernuclei on the example of j^^Ph have 
been calculated in relativistic DDRH theory. The initial state was represented by a ground state ^°^Pb core with 
an attached A single particle state, the final state as a neutron particle-hole excited ^°*Pb where the A hyperon 
occupies a lower lying single particle level and the neutron of the particle-hole pair being unbound. Hence the process 
corresponds to a decay of the initial single A configuration into a A-neutron-particle-hole configuration. 

Due to the fact that the hypcrnuclear Auger effect appears mainly in heavy and intermediate mass hypernuclei the 
spectral distribution of the emitted neutrons is extremely complex. For this reason a very detailed reconstruction of 
each event will have to be done in possible experimental measurements. By focusing on initial states with the A in 
the 1/^-orbit our calculations predict a clean doublet structure in the Auger neutron spectrum from which the 1/ 
spin-orbit splitting can be directly read off, assuming sufficient energy resolution. For a general A orbit no definite 
signature of the spin-orbit splitting will appear due to the huge amount of transitions and the additional broadening 
of peaks when deeply bound neutrons arc involved. The effect of a finite width of the neutron states, increasing with 
the distance from the Fermi-surface has not been considered in our calculations. For the least bound neutrons it will 
have almost no effect, but going deeper in binding energy the broadening of the states can no longer be neglected. 
As described in the previous paragraphs the neutron spectra in which the deeply bound neutrons are also involved 
are already fairly complex. An additional broadening of the neutron levels will thus wash out most of the spectral 
structure there. However, in the case of the l/'^-shell de-excitation only the valence neutrons close to the Fermi-level 
are involved so that the spin-orbit splitting signal is going to survive in realistic spectra. 

For the spectroscopy in the intermediate and heavy mass region the hypcrnuclear Auger effect provides a promising 
complementary tool to 7 spectroscopy, although special care must be taken on the choice of transitions and tagging 
the energy. For the ^°^Pb hypernucleus the 1/^ shell is such a well suited case. This technique might be the only 
way to achieve high resolution information on the hypcrnuclear fine structure in heavy nuclei since for 7 transitions 
it will be even harder to assign the detected photons to specific transitions. 
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APPENDIX A: MATRIX ELEMENTS 

For the evaluation of the matrix elements the transition form- factors, eq. (^, are evaluated in the spherical basis 
of Dirac wave functions: 

r is either the four by four unit matrix or 7^, depending whether the scalar or the vector potential is evaluated, q and 
f denotes the unit vectors in direction of q and r, respectively. The spinors are the usual total angular momentum 
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eigenstates of the Dirac equation with the generahzed spin-angle spherical harmonics Qjim{0, (j)) obtained by coupling 
spin and orbital angular momenta 24 . The numerical solution of the radial Dirac equation was discussed e.g. in 
[E0[. The orbital angular momenta I and I are determined by j and the parity tt: 



j + 1/2 for TT = 
j-1/2 for TT 

For the scalar vertex we get then 



I 



i-l/2 for TT = 
j + 1/2 for TT = 



(A2) 



P?2('?) = E(-)'^a;(9) {[/ 



dr gi{r)g2{r)jx{qr) 



dr fi{r)f2{r)j\{qr) 



dfl n* J YxM . r 



(A3) 



The radial matrix element is evaluated numerically. The angular integral can be performed analytically. By means of 
the Wigner-Eckhardt theorem pOt] it can be expressed by Clebsch-Gordan coefhcients and reduced matrix elements 



For the vector transition form factor we have to evaluate 

.(9) = E(-)'^a;('Z) {[/ drr^ g,{r)g2{r)n{qr) 



3 A j' 
-m fi m' 



2-' 



(A4) 



dr fi{r)f2{r)j\{qr) 
dr gi{r)f2{r)j\(qr) 
dr r^ .fi(r)g2{r)j\{qr) 



Slhmi J2t2m2 



(A5) 



The a are the usual Pauli matrices. The reduced matrix element for the angular matrix element is identical to 
equation (A4), for the spatial components it is given by 



{jlm\ Yx^aM \]'l'm') = ^ IN) {]lm\ \Yx^.aiM]iN \j'l'm') 



IN 



IN 

= E(-)'"" IN) 



IN 



J I f 

m N ml 



J I f 
m N ml 



2-' 



2-' 



I 1 J 

I I / }{i\\yx\\n 

A 1 / 



(A6) 



Note that the Pauli matrices ctm, M = 0, ±1, are used here in the spherical basis Explicit expessions for the 

reduced matrix elements are found e.g. in ref. pM. 
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